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Hydrodynamics with spontaneous symmetry breaking: application to relativistic 

heavy ion collisions 
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In this paper we apply hydrodynamics for systems with continuous broken symmetries to heavy 
ion collisions in the framework of (1+1) dimensional Bjorken model. The temperature profile with 
respect to proper time determined in that context exhibits no differences with the ideal fluid. On 
the contrary, it is shown that the profile obtained when Muller-Israel-Stewart second order theory of 
dissipation is included on top of standard hydrodynamics indicates a slower cooling of the system. 



I. INTRODUCTION 

Relativistic hydrodynamics is a very useful tool in the 
field of heavy ion collisions when one wants to reproduce 
the space-time evolution of the hot and dense matter cre- 
ated in such events. Upon the hypothesis of local equilib- 
rium, which can be questionable for such small and not 
long-lived systems, this approach directly gives in terms 
of usual quantities such as temperature, pressure..., an 
intuitive description of what happens during the whole 
evolution of the system until the freeze-out takes place. 
Moreover, it is usually assumed that the matter under 
consideration is almost perfect, so that standard hydro- 
dynamics can be applied directly. However, the hot and 
dense matter produced is mainly constituted of pions 
and, with the exception of no attention has been 
drawn to the fact that the breaking of chiral symmetry 
can affect the theory of hydrodynamics itself. One of the 
aims of this paper is thus to give in a simple approach 
a quantitative estimate of this modification. More pre- 
cisely, we will determine the temperature profile with re- 
spect to time in the framework of the (1+1) dimensional 
boost-invariant Bjorken model with and without the ef- 
fect of the spontaneous breaking of the symmetry In a 
second step we also introduce the dissipation process and 
compare the previous results with those obtained when 
Muller-Israel-Stewart second order theory of dissipation 
is included on top of standard hydrodynamics. 



central rapidity region of a heavy ion collision, we can 
neglect in a first step the baryonic density. Thus, the 
hydrodynamic equations taken from become: 



(1) 



iMif? ~ /^tt'VS&Et + / s 2 £^£t) = (2) 

where the energy-momentum tensor is given by 

T» v = {e+p)u^u l -pg» v + f 2 /4tT(dt i Ed u Y$ + d l, T,d^). 

In the above equations, ft (/ s ), which can be inter- 
preted as the temporal (spatial) pion decay constant, are 
actually not constants but have a non-trivial temperature 
dependence. Nevertheless, in order to have an estimation 
of the effects induced by the presence of the spontaneous 
breaking of the chiral symmetry, we will adopt two sim- 
ple approaches. The first one, for the sake of simplicity, 
is to consider that ft = fs = fix = 93 MeV. Of course, 
at finite temperature, this is not true and, for instance, 
for small temperatures, one knows that f 2 and f 2 differ 
from f 2 by a factor proportional to T 2 while the cor- 
rection for f 2 — f 2 is of order T 4 so that, to first order, 
f 2 = f 2 = / 2 (1— T 2 /(6/ 2 )). This second approach, more 
realistic, is based on this formula, extended to larger tem- 
peratures. 



III. THE BJORKEN PICTURE 



II. HYDRODYNAMICS EQUATIONS FOR 
CHIRAL FLUID 

Following pj, the hydrodynamics degrees of freedom 
in a chiral fluid with chiral SU(2) symmetry sponta- 
neously broken are the densities associated to conserved 
quantities, namely entropy density s, momentum density 
T° l , baryonic number density n, left and right-handed 
charge densities which can be written as ST/(2)-matrices 
Pl = p L Ti/2 and pr = p R Ti/2 and finally the variables 
associated to the Goldstone modes introduced in the the- 
ory via a SL/(2)-matrix £ = e ir ' 7r /^' r . If we focus on the 



Even if several important improvements concerning the 
geometric description of a heavy ion collision have been 
done since the pioneering paper of Bjorken |^|, we will 
use here this simplified version. The Bjorken picture of a 
collision is based on the hypothesis that initial conditions 
are invariant under Lorentz transformations throughout 
the central rapidity plateau. This implies that, for a 
one-dimensional flow in the z-direction, all the quantities 
depend actually not on the time t and the position z 
but only on the proper time r = \/t 2 — z 2 . The space- 
time evolution of the system can be then represented by 
hyperbola r = est. If the hydrodynamic velocity is 
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taken as w M = x^/t = (t/r, 0, 0, z/t), the hydrodynamic 
equations 



dyT^ = 



(3) 



with = (e+p)u^u u — pg^ v lead for an ultra relativistic 
gas to the following differential equations 0, : 



ds s 
dr t 
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whose solutions give the following scaling behavior 
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(4) 
(5) 
(6) 



s(t) = s(r )- 



(7) 
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T(r)=T(ro)(^) 1/3 (8) 
<r)=6(r )(^) 4/3 . (9) 



And, since p = e/3: 



p(r) =p(r ) (^J 



4/3 



(10) 



IV. RESULTS 



We first take into account the effect of the spontaneous 
breaking of chiral symmetry where no renormalization of 
the pion decay constant occurs, i.e. ft = fs = f-n = 93 
MeV. We thus obtain from the hydrodynamic equations 
(JTJ , J2J in the first non- vanishing order of the pion field : 



ds 



= o 
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9r t ) dr 



= 0. 



(11) 



(12) 



(13) 



(14) 



and from the explicit expression of the pressure p given 
in Q: 



(15) 



We see that there is no modification of the equation for 
the entropy density and for the temperature, that is the 
scaling behavior is exactly the same s(r) = s(t )t /t and 

1 /3 

T(r) = T(t ) (t /t) 1 . On the contrary, for the energy 
density and for the pressure, one can observe a modifica- 
tion due to a term proportional to dir 1 /dr which can be 
determined from the last equation (|14|) • The solution for 
the time derivative is d^/dr = A/t where A is a con- 
stant whose value is fixed by initial condition for r = tq. 
If we now report in the above equations, we can solve 
analytically the resulting system to obtain: 
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p=p 
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(16) 



(17) 



In order to fix the initial condition related to the value 
of A, let us now consider the density matrix p = exp~^ H . 
The part of the Hamiltonian relevant for chiral hydrody- 
namics degrees of freedom can be taken from [l|. An 
explicit calculation with Bjorken modelisation of the col- 
lision gives: 



p = exp[-/3 / d J xT 00 oc exp 
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Assuming isospin 
p cx exp {— ^ J 



d 3 : 



(18) 

symmetry, we then have: 

2^ 



^)"J- For t = To, we thus 
immediately see that p is Gauss distributed. Actually, if 
we call Vo the initial volume in the rest-frame of the two 
colliding nuclei and 7 cm the Lorentz factor needed to go 
from that frame to the laboratory frame, we have: 



/ 3/3 V ( A s 
p(r ) cx exp — — 

\ 2 7 cm \ r 0/ 



(19) 



We can easily estimate Vo by geometric arguments: Vo = 
itR 2 2tq where R is the radius of the colliding nuclei. If we 
now put typical numbers on all these constants: R ~ 5 
fm, r = 1 fm, 7 cm ~ 10 and T(r ) = 200 MeV, we 
can directly estimate the effect of the terms associated 
to the spontaneous symmetry breaking since < A 2 > is 
simply related to the Gaussian distribution: < A 2 >= 
T(T )T ~f cm /(6nR 2 ). Numerically, using these inputs, we 
obtain a correction of order ~1% with respect to the 
perfect case. 

We now introduce the explicit chiral perturbation T- 



dependence for f t and f s , i.e. f 2 



mi 



T 2 1 (6f 2 )). The evolution of the temperature is then gov- 
erned by the following equation : 



AcltT a - 



< A' 2 > 



r(l-T7(6/2)) 2 6/2 



s(t )t (20) 



where now < A' 2 >=< A 2 > (1 - T$/{6f 2 )). The nu- 
merical results are presented on figure ^ We can observe 
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that the correction, which induce a faster cooling of the 
system is of order 5%. 

Even if the approach can be considered too simplistic, 
it is not reasonably expected that an improved treatment 
can modify noticeably the order of magnitude. One can 
thus conclude that the effect due the spontaneous break- 
ing of chiral symmetry is quite negligible and without 
any experimental consequences. 

V. DISSIPATIVE EFFECTS 

In the previous part we have seen that it is possible 
to neglect the effect of the breaking of the symmetry on 
the hydrodynamics theory itself, that is we can use " nor- 
mal" hydrodynamics. Nevertheless an important ques- 
tion which has already been addressed in Q is the effect 
of dissipation on top of this model of perfect fluid. 

There are two non equivalent possibilities to take into 
account dissipation. First, one can use the usual formal- 
ism and just say that when an inhomogeneity appears in 
a system, the latter reacts and a flow proportional to the 
gradient develops itself in order to restore equilibrium. 
This is first order theory. However it is well-known that 
this approach presents instabilities and that it is neces- 
sary to consider a relaxation time for the flow. This goal 
is achieved through Miiller-Israel-Stewart second order 
theory. In the following we will explore both methods. 

In 0, the author showed that only the dissipative ef- 
fects related to the shear viscosity r\ play a role for the 
temperature profile in the (1+1) dimensional Bjorken 
model. The main physical reason is that bulk viscosity 
vanishes for an ultra relativistic gas. He took the shear 
viscosity from where the temperature dependence for 
X] is rj = irf£/(8T). Then, the differential equations to be 
solved for the energy density e = 3aT 4 are: 



where $ is related to the stress tensor (see for further 
details) and $ = for the ideal (perfect) fluid, <& = 
for the first order and ivS^ = — $ + |- for the second 

77 or 3 r 

order theory discussed before. In the last equation r n is 
the relaxation time associated to $. 

We now use (|21[l , using the same inputs as , namely 
massless pions and current algebra cross section. We can 
see on figure that the correction due to viscosity com- 
pared to the ideal case is totally negligible not only for 
the (non-physical but usual) first order but also for the 
second order. Our results are different from those pre- 
sented in the original paper . The difference lies in fact 
in a small computational mistake in . The main point 
is that, contrary to what is claimed in 0, there are, in 
this simple approach, no appreciable effects. If we would 
stop here, we could say, as in the first part dedicated to 
the breaking of chiral symmetry that a realistic descrip- 



tion of a pion gas can be achieved with the perfect fluid 
model. 

However, we stress that this conclusion only holds be- 
cause of the hypothesis used in |2|. The temperature 
dependences of the viscosity 77 and of the relaxation time 
7v are only valid for massless pions and when the cross- 
section for 7T— 7r interaction is given by current algebra. If 
we take instead the realistic case of massive pions and the 
experimental cross-section, these dependences are com- 
pletely different. For instance, r\ increases with temper- 
ature Q, |1| instead of decreasing. With this realistic In- 
dependence for the viscosity and the relaxation time for 
the pions taken from ||, we obtain the curves plotted 
on figure El We can now observe that the first and the 
second order are still close but they are now around 20% 
above the ideal fluid case. This means that dissipative 
effects imply a slower cooling of the system and should 
be incorporated in the simulations of heavy ion collisions. 



VI. CONCLUSION 

In the first part of this paper we have examined in a 
simple model an application of hydrodynamic equations 
with spontaneous breaking symmetry. In order to get an 
estimate of the magnitude of the effects, we have used 
the one dimensional Bjorken description of a heavy ion 
collision. We have shown that the model of the perfect 
fluid is almost correct : the effect of the breaking of the 
symmetry is to give a correction of about 1% when no 
T-dependence of the pion decay constant is taken into 
account and a few percents when we consider a crude 
T-dependence. In a second part we have considered an 
other deviation from the perfect fluid model, namely the 
dissipation. The dissipative effects, when taken into ac- 
count in a simple description 0] are not so important as 
claimed in this reference but represent actually a few per- 
cents. On the contrary, we have obtained sizable effects 
with a correct T-dependence for the shear viscosity and 
the shear flow relaxation time. For example, for an ini- 
tial temperature of 200 MeV, the temperature goes from 
about 110 MeV in the ideal case to 140 MeV with our 
dissipative effects for t/t ~ 6. The whole conclusion of 
this paper is that the temperature profile generated from 
the perfect fluid approach and used in many simulations 
is almost correct if one takes into account only the effects 
of symmetry breaking. On the contrary, if we consider 
dissipative effects the profile can be substantially modi- 
fied, that is we have a slower cooling of the pion gas. 
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FIG. 2: Temperature profile with respect to proper time for 
the perfect fluid (solid line), first order and second order the- 
ories. 



